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I. INTRODUCTION 


The Boussinesq equation, a cornerstone in the time series 


Abstract— This study develops a time series model for water waves 
using a weighted Taylor series approach to analyze water wave 
dynamics. The model incorporates the continuity equation, Euler's 
momentum conservation equation, and the Kinematic Free Surface 
Boundary Condition, all formulated through the weighted Taylor 
series. By integrating the modified continuity equation across the 
water depth, utilizing depth-averaged velocity concepts, we derive the 
water surface elevation equation. Similarly, applying the Euler 
momentum conservation principle to the water surface yields an 
equation for the horizontal water particle velocity, which is 
subsequently converted into an expression for the horizontal depth- 
averaged velocity. The equations for water surface elevation and 
horizontal water particle velocity are solved using numerical 
methods. The application of the numerical model results in the 
generation of four distinct wave profiles: sinusoidal, Stokes, cnoidal, 
and solitary wave profiles, classified according to Wilson's (1963) 
criteria. The emergence of each wave profile type is influenced by the 
specified input wave height. 


significantly minimizes, if not entirely eliminates, the 
truncation error typically associated with Taylor series 
approximations. 


water wave model, has been developed in various forms by 
numerous researchers, including Boussinesq (1871), 
Dingermans (1997), Hamm, Madsen, and Peregrine 
(1993), Johnson (1997), Kirby (2003), and both the 1967 
and 1972 works of Peregrine, among others. These 
versions encompass both the continuity and the water 
surface equations, as well as the water particle velocity 
equation. Each variant of the Boussinesq equation stems 
from the foundational principles of the continuity equation 
and Euler's momentum conservation equation, albeit with 
unique interpretations that lend each version its distinctive 
characteristics. 


In our research, we developed a time series model that 
refines the basic hydrodynamic equations—the continuity 
equation and Euler's momentum conservation equation— 
through a weighted Taylor series approach. This method 
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We derived the water surface elevation equation by 
integrating the continuity equation across the water depth, 
utilizing the concept of depth-averaged velocity. This 
integration employs the velocity potential equation from 
the solution to the Laplace equation for water particle 
velocity, ensuring precision in the integration process. 


Furthermore, the horizontal water particle velocity 
equation was formulated by applying Euler's momentum 
conservation principle at the water surface. This equation 
was then converted into a horizontal depth-averaged 
velocity equation, leveraging the horizontal water particle 
equation from the velocity potential solution of the Laplace 
equation, thereby achieving exact transformation. 
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The core equations of our study—the surface water 
elevation equation and the horizontal depth-averaged water 
particle velocity equation—were solved numerically. The 
spatial differential was tackled using the Finite Difference 
Method, while the time differential was addressed through 
the predictor-corrector method. Both methods are 
characterized by second-order Taylor series accuracy, 
contingent upon the selection of an appropriately small grid 
size and time step to ensure the Taylor series is truncated 
to second-order differentials only. Our research calculated 
the grid size and time step to maintain the Taylor series 
truncation to first-order, aligning with the accuracy 
requirements of the Finite Difference Method. 


Meticulous formulation of the basic hydrodynamic 
equations, governing equations, and numerical solutions, 
all executed with high precision will generate more reliable 
and accurate results. 


Il. WEIGHTED TAYLOR SERIES 
Weighted Taylor series is a Taylor series that is truncated 
to only a first-order series, where the contribution of 
higher-order terms is expressed with coefficients called 
weighting coefficients (Hutahaean (2023a)).. 


For functions with two variables f = f (x, t) 
f(x + 6x,t + ôt) 


d d 
= ft) + yet A+ yðr nn. (1) 


Yt,2 and y, are weighting coefficients. 
For function f = f(x,z,t) the weighted Taylor series is 


d 
f(x + 6x,z + 6z,t + ôt) = f(x,z,t) + vade L 


af af 
+yxôx at YzÔôZ apo (2) 


Yt3> Yx and y; are weighting coefficients. There is no 
difference between yy in f(x,t) and yyin f((x,z,t). The 
basic values of the weighting coefficient are, 

M42=23 W3=35%=135%,=1 


The foundational values in our model are obtained by 
eliminating higher-order even differential terms using the 
central difference method, a technique detailed by 
Hutahaean (2003a). The adjustments to weighting 
coefficients, addressing odd higher-order differential 
terms, are outlined in Table 1, where these coefficients are 
defined in relation to the optimization coefficient, €. This 
relationship indicates that an increase in e directly 
amplifies the impact of higher-order Taylor series terms on 
the analysis. Although the methodology for determining 
these coefficients is elaborated in Hutahaean (2023a), it is 
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important to emphasize that the accuracy of the coefficients 
listed in Table 1 significantly exceeds those previously 
reported by Hutahaean (2023a), reflecting advancements in 
the precision of incorporating higher-order differential 
terms to enhance the model's predictive accuracy and 
reliability in hydrodynamic evaluations. 

Table (1) Weighting coefficients. 


E Yt,2 Yt,3 Yx Yz 
0.010 | 1.9998 | 3.0049 | 0.9988 | 1.0115 
0.020 | 1.9992 | 3.0210 | 0.9951 1.0490 
0.022 | 1.9990 | 3.0257 | 0.9941 1.0602 
0.024 | 1.9988 | 3.0310 | 0.9930 | 1.0726 
0.026 | 1.9986 | 3.0368 | 0.9918 1.0866 
0.028 | 1.9984 | 3.0433 | 0.9904 | 1.1020 
0.030 | 1.9981 | 3.0503 | 0.9890 | 1.1189 
0.032 | 1.9979 | 3.0581 | 0.9875 1.1375 
0.034 | 1.9976 | 3.0664 | 0.9858 1.1577 
0.036 | 1.9973 | 3.0754 | 0.9841 1.1798 
0.038 | 1.9970 | 3.0851 | 0.9823 1.2037 
0.040 | 1.9966 | 3.0955 | 0.9803 1.2296 


€ is the optimization coefficient, where in this time series 
model the selection of the weighting coefficient is 
determined by the wave amplitude the greater the wave 
amplitude the greater the €, which is in the range of 
0.015 < £e < 0.035. 


II. DEPTH INTEGRATION AND 
TRANSFORMATION COEFFICIENT 
The formulation of the water surface equation is done by 
integrating the continuity equation with respect to water 
depth. The integration is carried out using the concept of 
depth average velocity (Dean (1991)) as follows. 


-im 
a S , (3) 


U = U(x,Zo,t) is the horizontal particle depth average 
velocity, Fig (1) 

Zo is the elevation of U towards still water level. 

u = u(x, z, t) is the horizontal water particle velocity 

Pu is integration coefficient 

D is total water depth, D = h +n 

h is water depth towards still water level 

7 is surface water elevation towards still water level 
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Ai water level 


U = u(x, Zo t) 


Sea bed 


Fig (1). Depth average velocity concept 
a. Integration coefficient p, 
In this study, a method for calculating the integration 
coefficient f, was developed using the velocity potential 
equation and by establishing that the depth average 
velocity is U = u(x,Zọ,t), where the value of Zo is 
determined (see Fig (1)). 


Velocity potential solution of Laplace equation, at the 
characteristic point where cos kx = sinkx (Hutahaean 
(2023b)) is, 

o(x,z,t) = 2G cos kx coshk(h + z) sin øt 

(x, Z, t) is potential velocity 

G is wave constants 

k is wave number 

h is water depth towards still water level 

o is angular frequency 

The horizontal velocity equation u(x,z,t) is obtained 
using the velocity potential property: 


u(x, z,t) = -22 
u(x,z,t) = 2Gk cos kx cosh k(h + z) sin ot (4) 


It has been defined that the horizontal depth average 
velocity is the velocity at z = Zọ (Fig (1)), where Zp is 
determined, then 

U(x, Zo, t) = 2Gk cos kx cosh k(h + Zo) sin ot 


Employing (4) and (5), 
u(x,z,t) _ coshk(h+z) 
U(x,Zo,t) cosh k(h+zo) 
Equation (3) is written into the equation for the integration 
coefficient p, 


1 m 
Bu = Gp Pia 
1 n 
= —— h k(h d 
Bu Deoshktez) | °° ee ee 


Integration is completed, 
_ sinh k(h + 7) 
4" kD cosh k(h + zo) 


As 3 
=> A is wave amplitude, therefore 


pane 
TERZ 
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sinh kD 
Bu = kD coshk(h + Zo) 
Hutahaean (2024) found that kD = Om where tanh kD = 
tanh 6m ~ 1, @ is the deep water coefficient, in this time 
series model research used 0 = 1.80. Determined Zọ = 
—th, where 0 < € < 1.. In deep water, h ~ h + Z, 


sinh 07 


Bu=Geesh(onG=D) e (7) 


b. Transformation Coefficient a, 

In the process of integrating the continuity equation with 
respect to depth, there will be a horizontal surface water 
particle velocity u,. This velocity must be transformed into 
the horizontal depth average velocity U(x, Zo, t). 

Using (4) and (5), 

un cosh k (h + 5) 

U coshk(h+ Zo) 

Un cosh 67 

uU cosh(@(1 — é)) 

Transformation coefficient a, is 


cosh 61 
Xun = cosh(62(1-)) -++(8) 
Un = AynU ..(9) 


In the formulation of the equation for horizontal depth 


average velocity, the variable u,u, must be transformed 


n 
into U. 


Using (4), 

Uy = 2Gk cos kx cosh k(h + n) sin ot 

UyUy is mathematically expressed as 

Unun = ( 2Gk cos kx cosh k(h + n) sin ot)? 

However, since u,U, moves simultaneously with u, 


where the distribution of u,,u, according to space and time 


n 
is the same as u,,, thus, 

UyUn = (2Gk)* cos kx cosh k(h + n) sin ot 

UU = (2Gk)? cos kx cosh k(h + Zo) sin ot 

Therefore, the relation between u,u, and UU is 

Ugly = AUU ee (10) 
As an illustration of the values of the integration 
coefficient f, and the transformation coefficient «uy can 
be seen in Table (2). 


Table (2) The values of By and aun 
Š Pu Gun 
0.10 0.324 1.679 
0.15 0.420 2.176 
0.20 0.544 2.819 
0.25 0.705 3.653 
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0.30 0.913 4.732 
0.35 1.182 6.130 
0.4 1.531 7.937 


Table (2) shows that larger € is followed by deeper Zo, the 
larger the values of the integration coefficient f, and the 
transformation coefficient @,,,. In this study, a value of 
Pu = 1 was used, and using the Newton-Rhapson method, 
č = 0.317612 was obtained with a transformation 
coefficient value of @,, = 5.183954. 


IV. WATER SURFACE ELEVATION 
EQUATION 
4.1.Weighted continuity equation. 
The continuity equation is formulated using the weighted 
Taylor series and by applying the principle of mass 


conservation, 


d d 
Vests =O ast (11) 


This equation is formulated under the condition that the 
horizontal particle velocity only changes along the 
horizontal axis, and similarly, the vertical particle velocity 
only changes along the vertical axis. 


4.2. Integration of continuity equation 

To obtain the water surface elevation equation, the 
continuity equation is integrated with respect to water 
depth, where the integration of the second term can be 
directly completed by ignoring the bottom vertical water 
particle velocity. 


n du 
[rade trem =0 


The vertical surface water particle velocity w, is 
substituted with the Kinematic Free Surface Boundary 
Condition, where the weighted KFSBC (Hutahaean 
(2023a)) is as follows. 

dn dn 
Wy = Yt2 77 dt + Yxn ay 


Substituted into the integration of the continuity equation, 


n du dn dn 
| ra” + Yz (rz art VxUn =) =0 


The integration of the first term is solved by Leibniz 
integration (Protter, Murray, ie & Charles, 1985), 


Ba B 
i ad a eT (12) 
1 du dn dh 
| zdz -Ef u dz — Uy — Unz 
-h dx dx dx 


The integration of the first term on the right side of the 
equation is completed using the concept of depth average 
velocity, and the horizontal bottom water particle velocity 
is ignored. 
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1 du df„UD dn 
z l=- 
-n ax dx dx 
Substituted into the continuity equation integration, the 
surface velocity is transformed into depth average velocity, 
and the equation is written as the water surface equation, 


dn _ dUD dn 


VzVt,2 Fe dt —VxBu az + Vx (1 — Vz) Qy,U ax 


(13) is the water surface elevation equation that is used to 
calculate the water surface elevation. 


V. EULER’S MOMENTUM CONSERVATION 
EQUATION 

By using the weighted Taylor series and the same fluid 
flow conditions as those used in the formulation of the 
continuity equation, where horizontal water particle 
velocity only changes along the horizontal axis and vertical 
water particle velocity only changes along the vertical axis, 
Euler’s momentum conservation equation in the horizontal 
and vertical directions, respectively are presented as 


follows. 

du du _ z 1 dp 

at ea oie (14) 
dw dw _ _ 1dp = 

a aS, nae g (15) 


Equation (15) is written as an equation for pressure p and 
integrated over water depth and the surface dynamic 
boundary condition ss = 0 is performed 

p "dw 

-= rea | — dz +2 5 (ww — ww) +g -2z) 

p z dt 

Differentiated with respect to the x-horizontal axis, 


1dp d dw yz d 
Pong ae Zt ga aw) 
dn 
KET 


Substituted to (14) 
du Yx duu _ d F dw 


Yesae 2 ax ax] ae 
-2E (w W,W, — ww) — g? 


2 dx 


This equation is worked out on the surface i.e. on Z = 7, 


du, Yx duyUy dn 
Yt,3 = -g 
dt 2 dx dx 
The surface velocity is transformed into depth average 
ee 
Yx@un AUU a 
Viton x ae as =-g7 PER (16) 
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VI. VERTICAL WATER PARTICLE VELOCITY 
EQUATION 
Both the water surface elevation equation and the 
horizontal water particle velocity do not include the 
variable of vertical water particle velocity. However, this 
velocity can be calculated in a simple manner, using the 


Kinematic Free Surface Boundary Condition: 
_. dn r dn 
Wn = Yez Ge Ma ae 
This surface velocity is transformed into depth average 


velocity. 


EPE an an 
w= (Vaz + aun U 2) Ds 
Where, 
sinh 07 
awy = shen) ...(18) 


Where @ and € are presented in Section 3. 


VIL. NUMERICAL SOLVING 
Equations (13) and (16) are solved using numerical 
methods. The spatial differential is solved with the Finite 
Difference Method, and the time differential is solved 
using the predictor-corrector method. 


7.1. Calculation of time step dt and grid size 6x 

Solving differential equations with numerical methods 
requires a specific time step and grid size to ensure good 
calculation results. The following section presents the 
equations for calculating these time steps and grid sizes. In 
principle, these equations are formulated based on the 
truncation of the Taylor series, which is the interval size 
where the Taylor series can be truncated to only a first- 
order series. The formulation of the equations can be seen 
in Hutahaean (2023a). 


a. Time step dt calculation 

Time step ôt was measured using the following equation. 

“ne? — ne, te =0 ....(19) 
€ is a small number that determines the size of the time step 
ôt for example, € = 0.01. The smaller the value of e, the 
smaller €, and the smaller dt, where ôt = €T, T, with T 
being the wave period. In that equation, there are two 
values of £+; the smallest value is used. 


b. Grid size 6x measurement 

The calculation of grid size 6x requires the input time step 
ôt obtained from (19). The equation to calculate ôx is, 
Cot Cex + C2 E = .....(20) 


Co = 2Te,E — 2r? e? 
Cy = 2me + 4772, 


Cp = —2n? 
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There are two values of £, the largest is used. Grid size 
ôx = €,L,, where L is the wave length. 

The formulations (19) and (20), with higher accuracy, can 
be seen in Hutahaean (2023a). However, the dominant 
factor determining accuracy is the choice of the value of e. 
In principle, the smaller the value of £. the better. In Table 
(3), the calculation results for ôt and ôx for waves with a 
wave period T = 8.0 sec, and a wave length L = 16.0 m 
are presented. 

Table (3). Time step ôt and grid size 5x outcomes 


ôt ôx 


E Et Ex (sec.) (m) 
0.02 | 0.0065 | 0.0192 | 0.0516 | 0.3079 
0.025 | 0.0081 | 0.0241 0.0648 | 0.3856 
0.03 | 0.0097 0.029 0.078 | 0.4637 
0.035 | 0.0114 | 0.0339 | 0.0913 0.542 
0.04 | 0.0131 | 0.0388 | 0.1047 | 0.6207 
0.045 | 0.0148 | 0.0437 | 0.1183 | 0.6997 
0.05 | 0.0165 | 0.0487 | 0.1319 | 0.7791 


As seen in Table (3), the larger the €, the larger also the ôt 
and ôx. To check the compatibility between the time-step 
dt and the grid size 6x, the Courant criterion (1928) is used 


as follows La 3C 
ôt 


C adalah wave celerity dimana 
o 


L 
C= =F m/sec 


For example, for waves with a wave period T = 8 sec., and 


with a wave length L = 16.0 m, C = 2.0 m/sec. The 
ôx 


lar for ôt and 6x in Table (3) are calculated, 


with the calculation results presented in Table (4). 
Table (4). Checking the Courant criteria 


values of 


E ôt ôx T or 
(sec.) (m) (m/sec) C 

0.02 | 0.0516 | 0.3079 | 5.9639 2.982 
0.025 | 0.0648 | 0.3856 | 5.9547 | 2.9774 
0.03 0.078 | 0.4637 | 5.9454 | 2.9727 
0.035 | 0.0913 | 0.542 5.9361 2.968 
0.04 | 0.1047 | 0.6207 | 5.9267 | 2.9633 
0.045 | 0.1183 | 0.6997 | 5.9172 | 2.9586 
0.05 | 0.1319 | 0.7791 | 5.9076 | 2.9538 
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In Table (4), it is found that the values of dt and 6x meet 
the Courant criterion, for all values of £, which is ~ 3. 


Modeling in non-uniform water depths requires a non- 
uniform grid size as well. This is because different water 
depths have different wavelengths. In a non-uniform grid 
size, there is no difference or change in the Courant 
number. 

ox E&L & 

be aT ey 

The Courant number in this case is = This value is 
constant for a given accuracy level £, unaffected by 
wavelength or wave period. As an example, calculations 
for a sloping bottom with nonuniform wavelength are 
presented in Table (5), where a wave period of 8 sec is 
used, and the grid calculation uses ¢ = 0.02. 


In Table (5), ôt is constant, while ôt is not, but the resulting 
Courant number, 
is constant. 
Eat _ ex 
C Et 


Table (5) Grid-size of non uniform water depth 


L ôt 5x Oe) Sr/a 
(m) (sec.) (m) (m/sec) C 

16 0.05163 | 0.3079 5.9639 | 2.98195 
14 0.05163 | 0.26941 | 5.21842 | 2.98195 
12 0.05163 | 0.23093 | 4.47293 | 2.98195 
10 0.05163 | 0.19244 | 3.72744 | 2.98195 
8 0.05163 | 0.15395 | 2.98195 | 2.98195 
6 0.05163 | 0.11546 | 2.23646 | 2.98195 
4 0.05163 | 0.07698 | 1.49098 | 2.98195 
2 0.05163 | 0.03849 | 0.74549 | 2.98195 


The Finite Difference Method is formulated using a 
second-order Taylor series, where the Taylor series is 
truncated to second order only. Therefore, calculations 
with this method require small time steps and grid sizes. 
The time step and grid size calculated for the truncation of 
the first-order Taylor series definitely meet the 
requirements for the truncation of the second-order Taylor 
series. 


c. Wave length measurement 


For calculating grid size, information about the wave 
length is required. The following section discusses the 
method for calculating wave length. 
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The deep-water wave number is calculated using the 
equation, 


gAo k2 — g tanh(@7) 
29 Ve 


g gravitational acceleration 


ko + Ve2Ve30° = 0 ssl) 


Ao deep water wave amplitude 


o= =, T adalah wave period. 


Deep water depth, 
=e ou 
ho = ae vaa) 
Wave number of one water depth hy < ho, 
koh 
ky, = oho 
hy 
Where wave length 
a 21 
x ky 


To conform to the wavelength model, in (21) the weighting 
coefficients are used: yY¢2 = 1.8, Y¢,2 = 2.6, Yx = 1.0 dan 
Yz = 1.0. 


7.2.Predictor-corrector method 


When performing calculations with the finite difference 
method, a grid-points is established within the domain, 
using the grid size as previously discussed. At a given point 
i on the grid, the water surface elevation equation and 
horizontal water particle velocity can be rewritten as 
follows: 


dn; 
= F(t 
a = ie) 
dU; = G,(t) 
dat! 
a. Predictor 


The steps of predictor method requires the use of central 
difference method, 

Niprea = + 26tFi 

After all Wi tesd values are obtained, the predictor step for 
water particle velocity is performed, 

Upsets = UE + 26tG} 


t+dt dan yt+ôt 


ipred iprea are obtained, the 


After all the values of 7 


corrector step is performed. 


b. Corrector 

Corrector is done by using numerical integration of 
Newton-Cote with 3 integration points (Arden, Bruce W. 
and Astill Kenneth N. ,1970), as follows. 


t+ôt 1 4 1 
f(t) dt = ôt (srt + -f? +4 sft) 
t-dt 3 3 3 


The three-point numerical integration is formulated using 
a three-point Lagrange polynomial that has an accuracy 
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equivalent to a second-order Taylor series. The time step 
calculated by (19) satisfies the time step requirement where 
the Taylor series can be truncated to second order. 

Where nf*2!, dan Uf*5", dihitung FE*5t 


ipred i,pred ipred 
t+ôt _ „t-ôt 1 t—ôt 4 t 1 t+ét 
Nicor = Ni + ôt ahi T ahi ae 3 Fiprea 


After all the He values are obtained, the corrector step 


for the water particle velocity is performed, 


1 4 1 
ma =o + 6¢(zaf% + Gf + a) 
After all values of fio! dan Uf%S*) are obtained, 


convergence is checked, |nji5t — ništa] <u, for all 
points, being a small number. If the convergence 
condition is not met, the result of the corrector step 
becomes the predictor result, Sol, = niet and 


ytret = ytret 


ipred icor> » and the corrector step is repeated. 


Generally, convergence is achieved in 3 replications 


VII. MODEL OUTCOME 

An assessment or evaluation of the modelled results 
was carried out by comparing the modelled wave profiles 
with the wave profiles from Wilson (1963). 
8.1. Wave profile based on Wilson's criteria (1963). 
The grouping of water wave profiles according to Wilson 
(1963) is presented in Table (6), with the profile sketched 
in Fig (2). 


trough 


Fig (2). Wave profile based on Wilson (1963) criteria. 


Table (6) Water wave profile criteria (Wilson (1963))) 


Wave type Nmax 

H 
Airy/sinusoidal waves < 0.505 
Stoke’s waves 0.505 — 0.635 
Cnoidal waves 0.635 — 1 
Solitary waves =1 


a. Airy’s Profile 
The Airy wave profile is marked by its symmetrical nature, 
with the wave crest and trough being mirror images of each 
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other, allowing for the wave height (H) to be approximately 
twice the wave amplitude (A), or H ~ 2A. This sinusoidal 
wave theory, also referred to as Airy wave theory or small 
amplitude wave theory, was pioneered by Airy in 1841. 
True to its designation, the theory is specifically tailored 
for waves of very small amplitudes, providing a 
foundational understanding of wave dynamics within this 
limited scope. 

b. Stoke’s Profile 

Stokes wave theory, introduced by George Stokes in 1847, 
represents an advancement beyond the sinusoidal or Airy 
wave theory, earning it the alternative name of Finite 
Amplitude Wave Theory. This theory addresses some of 
the limitations of the Airy wave theory by accounting for 
waves with larger amplitudes. The distinguishing feature 
of the Stoke’s wave profile is its asymmetry: the wave 
crests are sharper or steeper compared to the more gently 
sloped wave troughs, and it's possible for the wave profile 
to exhibit two distinct crests. 

c. Cnoidal wave profile 

The cnoidal wave theory was initially formulated by 
Korteweg and de Vries in 1895, upon observing the wave 
profile within a canal, thereby confirming its presence in 
natural settings. This particular wave profile is 
characterized by a significant disparity between the crest 
and the trough. Predominantly, the wave's surface is 
elevated above the level of still water, presenting a steep 
profile, whereas only a minor segment falls below, 
exhibiting a more gradual profile. 

d. Solitary wave profile 

In the case of a solitary wave profile, the entire surface of 
the wave remains above the level of still water. This type 
of wave was first identified by John Scott Russell in 1844 
after he observed the phenomenon in a laboratory setting. 
The analytical foundation for understanding solitary waves 
was later established by Joseph Boussinesq in 1871. Hence, 
like the cnoidal wave, the solitary wave profile is also 
observed in natural environments. 


8.2. Modeled wave profile in deep water. 

The model operated at a constant water depth of h = 30 m. 
At x = 0, a sinusoidal wave is introduced using the 
equation n(0,t) = A sin øt, where the wave period T = 
8 sec. This input is applied for the duration of one wave 
period only. The wave's amplitude is selected to match the 
desired wave profile. Calculations for determining the deep 
water depth hy, using equations (21) and (22), categorize 
the 30 m water depth as deep water, as indicated in Table 
(7). 

In the wavelength calculation results presented in Table 
(7), an inverse relationship is observed between wave 
amplitude and wavelength: as the wave amplitude 
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increases, the wavelength decreases. This trend contrasts 
with the findings. 


from the numerical model, where an increase in wave 
amplitude results in a slightly longer wavelength. It is 
important to note, however, that the water depth of 30 m 
exceeds hg and thus falls into the category of deep water. 


Table (7). Deep water depth hy measurement 


For the execution of this numerical model, a deep water 
coefficient 0 = 1.8. is employed. Additionally, a range for 
the weighting coefficient € = 0.015 — 0.035 is used, 
where a larger wave amplitude corresponds to a higher 
value of € the weighting coefficients are detailed in Table 
(1). The model runs over eight wave periods, equivalent to 
80 seconds. During the initial two wave periods, the wave 
profile undergoes an evolution, rendering it unstable. From 
the third wave period onward, the wave profile achieves 
stability. 


a. Sinusoidal profile 
With input wave amplitude A = 0.02 m, oe = — = 0.5 


is obtained (Fig(3)). This condition corresponds to 


Wilson's criterion for sinusoidal wave profile. 


0.03 
0.02 = 
2 
0.01 © 
= H = 0.04m Il 
> 0 mne RAA = 
` g 
= 0.01 = 
-0.02 
-0.03 
110 130 150 170 190 
x (m) 
eta (m) ——U (m/sec) W (m/sec) 


Fig (3). Sinusoidal wave profile 


b. Stoke’s wave Profile 
Stoke's profile, obtained using A = 0.5 m. The modelled 


results can be seen in Fig (4). Tmax = ae = 0.605 is 


obtained, this condition is in accordance with Wilson's 
criteria for Stoke's wave profile. 
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Fig (4). Stoke’s profile 
c. Cnoidal wave Profile 


The model was executed with an input of A = 1.30 m, or 
H = 2.60 m. This is the largest wave height for a wave 


2.1 f : 
Imar T 0.808 was obtained, in 


period of 8 sec. 
accordance with Wilson's criteria, for a Cnoidal profile. 


The image of this cnoidal profile is presented in Fig (5). 
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> se H =(2.6m 
> 0.5 
= o AA 
as 
90 110 130 150 170 
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eta (m) ——U (m/sec) W (m/sec) 


Fig (5). Cnoidal profile 


In this section, it was found that the model execution in 
deep water did not obtain a solitary wave profile, even 
though the maximum wave amplitude was used. 


8.3. The Effect of water depth on wave profile. 

In this section, the effect of water depth on wave profiles 
is explored using a model that simulates a sinusoidal wave 
input at x = 0.0 expressed as 7(0,t) = A sin øt, which is 
applied for the duration of a single wave period of 12.0 
seconds. The model uses a wave amplitude of 2.0 meters, 
equivalent to a wave height H = 4.0 m. For this scenario, 
the deep water depth hy is calculated to be 41.0 meters. The 
experiment is conducted at water depths of 41.0 meters, 20 
meters, and 8.0 meters to observe the resulting wave 
profiles. At the depth of 41.0 meters, matching the deep 
water depth ho, the model produces a Stoke's wave profile, 
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depicted in Fig (6). When executed at a water depth of 20 
meters, the model forms a cnoidal wave profile, as shown 
in Fig (7) and in a water depth 8.0 m, solitary wave profile 
was obtained (Fig.(8)). 
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Fig (8).Solitairy profile Nae = H, Amas =1 
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In the observed solitary wave profile, an increase in wave 
height is noted, from the initial 4.0 meters (input) to 4.40 
meters. This increase is believed to be a result of the wave 
setup phenomenon, which is characteristic of shallow 
water effects. According to Wilson's criterion for solitary 


Nmax _ 1 

H : 
However, the wave trough is not considered in this analysis 
due to its minimal size and the presence of numerous other 


wave crests. 


waves, the ratio of the maximum wave height 


8.4. Model execution at the sloping bottom. 

In a subsequent experiment, the model is applied to a canal 
featuring an upstream water depth of 15.0 meters and a 
downstream depth of 5.0 meters, with the channel 
extending over a length of 200.0 meters. The input for this 
setup is a sinusoidal wave characterized by a wave period 
of 8.0 seconds and a wave amplitude of 0.50 meters. The 
outcomes of this experiment are depicted in Fig (9), 
illustrating that at a water depth of 6.0 meters, the wave 
undergoes shoaling, which results in the wave height 
increasing to 1.40 meters and the wave adopting a solitary 
profile. 


The results from executing the model reveal its capability 
to simulate various wave profiles, including sinusoidal, 
Stoke's, cnoidal, and solitary profiles. Importantly, the 
formation of these profiles is influenced not only by the 
wave amplitude but also significantly by water depth, 
particularly in the case of the solitary profile. An amplitude 
that results in a Stoke's profile in deep water conditions 
may lead to a cnoidal profile in shallower waters, and 
similarly, the same amplitude can generate a solitary 
profile in even shallower depths. This underscores the 
critical role that both water depth and wavelength play in 
determining the specific wave profile that emerges. 


2 


1.5 


145 155 165 175 185 


x (m) 


eta (m) =———= U (m/sec) 


Fig (9). Solitary wave profil due to shoaling. 


The occurrence of a solitary profile in shallow water 
suggests that the model's wavelength is sufficiently short 
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to facilitate the formation of such a profile. The weighting 
coefficient emerges as a crucial parameter in shaping the 
wavelength, indicating that this coefficient, as part of the 
weighted Taylor series, significantly influences the 
outcomes of the model. This observation highlights the 
interconnectedness of wave amplitude, water depth, 
wavelength, and model parameters like the weighting 
coefficient in the complex dynamics of wave profile 
formation. 


IX. | CONCLUSION 


The study reveals that the time series model developed 
can generate wave profiles aligning with previous research, 
notably including the cnoidal and solitary wave profiles. 
This demonstrates the model's capability to effectively 
simulate water waves. 


The solitary wave profile emerges at short wavelengths, 
where the model's wavelength determination hinges on the 
weighting coefficients within the weighted Taylor series. 
This suggests the selected weighting coefficients are well- 
suited for the task. 


Moreover, the model's performance is influenced by the 
chosen time step and grid size. These parameters are 
aligned with the foundational hydrodynamics equations' 
assumptions, which presuppose truncation of the Taylor 
series to first order. Additionally, the time step and grid 
size adhere to the Finite Difference Method's requirements, 
where the Taylor series truncation is limited to the second 
order, ensuring a comparable level of numerical integration 
accuracy. The coherence between the time step and grid 
size also meets Courant's criteria, further validating the 
model's configuration. 


In summary, the employment of a weighted Taylor series 
to formulate the time series model, coupled with the 
strategic selection of time step and grid size for numerical 
solution, leads to robust model outcomes. 
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